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Abstract 

First results from RHIC on charged multiplicities, evolution of multiplicities with central- 
ity, particle ratios and transverse momentum distributions in central and minimum bias col- 
lisions, are analyzed in a string model which includes hard collisions, collectivity in the initial 
state considered as string fusion, and rescattering of the produced secondaries. Multiplici- 
ties and their evolution with centrality are successfully reproduced. Transverse momentum 
distributions in the model show a larger pT-ta.il than experimental data, disagreement which 
grows with increasing centrality. Discrepancies with particle ratios appear and are examined 
comparing with previous features of the model at SPS. 
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With the first colhsions at the Relativistic Heavy Ion CoUider (RHIC) at BNL in 
June 2000, the study of nuclear colhsions has entered the truly ultrarelativistic domain. 
While there exist predictions from many models |T[], now experiments have presented 
results 0, 0, |, |, ^, 0, |, ^ 0, [1^, [1^ on several aspects of data, most of them corre- 
sponding to AuAu collisions at 130 GeV per nucleon in the center of mass. So it comes 
the time to examine the ability of models for ultrarelativistic heavy ion collisions, fitted 
to describe nuclear data at the much lower energies of the Super Proton Synchrotron 
(SPS) at CERN and nucleon data in the range of energies going from SPS to TeVatron 
at FNAL, to describe the new situation, and whether the evidences of Quark Gluon 
Plasma (QGP) already obtained at SPS are verified or not |jT^. The aim of this letter 
is to compare the results of the String Fusion Model (SFM) |jl5|, |16|] with some of the 
first RHIC data. Other comparisons can be found in [|T^, |TB|Q After a very brief model 
description, charged multiplicities at midpseudorapidity in central collisions, evolution 
of charged multiplicities at midpseudorapidity with centrality, transverse momentum 
distributions of charged particles at different centralities and ratios of different parti- 
cles will be compared with available data coming from the experiments. Finally some 
conclusions will be summarized. 



An exhaustive description of the model can be found in Its main features are 
the following: Elementary inelastic collisions (binary nucleon-nucleon collisions) are 
considered as collisions between partons from nucleons of the projectile and the target, 
distributed in the transverse plane of the global collision. Some of these elementary 
collisions are taken as hard ones, and proceed as gluon-gluon — > gluon-gluon through 
PYTHIA [|19l with GRV 94 LO parton density functions (pdf 's) p| and EKS98 mod- 



ification of pdf's inside nuclei [21], with subsequent radiation and fragmentation per- 



formed by ARIADNE and JETSET ||T9[. Those collisions not being considered 
hard produce soft strings in pairs. These strings are allowed to fuse if their parent par- 
tons are close enough in impact parameter [|15[; as the number of strings increases with 
increasing energy, atomic number and centrality, this mechanism accordingly grows 
in importance. Fragmentation of soft strings is performed using the tunneling mech- 
anism for mass and transverse momentum distributions, while longitudinal momenta 
are simulated by an invariant area law. The main consequences of string fusion are a 

^In | [l8| a model which, hke ours, contains multipomeron exchange, a hard component and rescat- 
tering of secondaries, but no string fusion, is shown to be able to reproduce the experimental data 
on elliptic flow. 
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reduction of multiplicities in the central rapidity region and an increase in heavy par- 
ticle production. The produced particles are allowed to rescatter (between themselves 
and with spectators nucleons) using a very naive model with no proper space-time evo- 
lution, whose consequences are a small multiplicity reduction, an increase in strange 
and multistrange baryons and nucleon annihilation. Some comments are in order at 
this point: First, partons which generate both soft and hard strings can be valence 
quarks and diquarks, and sea quarks and antiquarks, so the number of soft strings is 
not simply proportional to the number of wounded nucleons but has some proportion- 
ality, increasing with increasing energy, centrality and nuclear size, on the number of 
binary nucleon-nucleon collisions|^. Besides, only fusion of two strings in considered in 
the actual version of the model, and hard strings are not fused. Finally, the rescattering 
model is simplistic and has been included just to estimate the effects that such kind of 
physics could have and to tune the parameters of the model as an initial condition for 
a more sophisticated evolution; thus, results depending strongly on it should be taken 
with great caution. All these aspects will be commented more extensively when the 
comparison with experimental data is performed. 

In Fig. |I| results of the model (unless otherwise stated, results of the model corre- 
spond to its default version with the mentioned pdf 's and string fusion and rescattering, 
see |]16[) for the pseudorapidity distribution of charged particles in central collisions at 
SPS and RHIC are compared with experimental data. For central AuAu collisions at 
130 and 200 GeV per nucleon in the center of mass, the model successfully reproduce 
the data (the ratio of multiplicities at 200 and 130 GeV is 1.08 in the model, slightly 



smaller than the experimental value 1.14±0.05 measured by PHOBOS [|T3|); while at 
56 GeV it overestimates the PHOBOS results 0. Nevertheless, the situation at these 
energies is not clear: WA98 results [H at SPS lie above the PHOBOS data at 56 GeV, 
and far above NA49 data |^ (as extracted in [0]) at SPS; NA49 results on multiplic- 
ities in central PbPb collisions at SPS are in agreement with those from WA97 ^ . 



So it is difficult to conclude anything definitive on the evolution from SPS to RHIC, 
of multiplicities with increasing energy in the model. 



^Usually the soft contribution is taken as proportional to the number of wounded nucleons, while 
the contribution proportional to the number of binary nucleon-nucleon collisions is considered hard. 
Let us stress that this is a misleading (model dependent) statement: some proportionality with the 
number of binary nucleon-nucleon collisions is demanded by a basic requirement of the theory as 
unitarity, and has nothing to do with the soft or hard origin of these binary nucleon-nucleon collisions. 
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Recently it has been proposed that the evolution of multiplicities with centrality 
can be used as a tool to discriminate among several models for multiparticle production 
in high-energy nuclear collisions. In this way, models which consider saturation of 
either the number of partons in the wave function of the projectile and target or in the 
number of partons produced in the collision show a constant or slightly decreasing 
behavior of the multiplicity per participant (wounded) nucleon with increasing number 
of participants^. On the other hand, models which consider some proportionality with 
the number of binary nucleon-nucleon collisions based on the AGK cancellation 



being this proportionality already present in the soft component |T6|, |31|, ^ or 
only in the hard component show a behavior, with the multiplicity per partici- 
pant increasing with increasing number of participants, qualitatively or quantitatively 
compatible with data. The results of our model for the 75 % more central collisions at 
SPS and RHIC are shown in Fig. |^ and compared with experimental data. It can be 
seen that the model underestimates WA98 data at SPS, while it overestimates those 
from NA49, as could be expected from the discussion about Fig. 0, but the qualitative 
behavior seems correct. At RHIC the agreement with data is quite satisfactory. It can 
be seen that the inclusion of rescattering results in a slight decrease of multiplicities, 
while the influence of string fusion is relatively small at SPS but very important at 
RHIC and crucial for the agreement with experimental data. In our model it is this 



latter mechanism the one which plays the role of shadowing corrections in [pT] , p^ , pi 



parton saturation in [^, ^ or string percolation in [^. Concerning the limita- 
tion of fusion of just soft strings in groups of two, let us point out that it seems to be 
compensated at RHIC with the choice of the fusion strength, while the non-inclusion 
of fusion of hard strings is unimportant, as they amount for just 1 % of the total num- 
ber of elementary inelastic collisions. This is no longer the case for the future Large 
Hadron Collider (LHC) at CERN, situation for which we present the results of the 
model in Fig. |] (results with rescattering are not presented because this mechanism 
is too CPU-time consuming at LHC energies for large nuclei): Here, the fusion of just 
two strings has reached its limit, so multiplicities are not so strongly damped as at 
RHIC, and fusion of more than two strings (and of hard strings, which now amount for 
32 % of the total number of elementary inelastic collisions), or even a phase transition 
like percolation ||3^, have to be introduced in the model. 

Other proposals which include saturation show an increasing behavior compatible with data. 
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Let us now turn to the transverse momentum spectrum. Preliminary measurements 
0, P] show that the spectrum in AuAu colhsions at 130 GeV per nucleon in the center of 
mass falls with increasing pj- faster than predictions from models which reproduce 
the pT-distributions in pp collisions at 200 GeV in the center of mass; this discrepancy 
grows with increasing centrality. A possible explanation is jet quenching i.e. the 
energy loss of high energy partons in a hot medium containing free color charges. So, 
there has been a great debate on the explanation of the absence of jet quenching at 



SPS and its presence at RHIC |3^, and its interpretation as a QGP signature. In our 
model we find quite the same feature as in , see Fig. ^, namely an excess of particles 
with high pt compared with experimental data, excess which becomes less pronounced 
when going from central to minimum bias collisions. Our model correctly reproduces 
multiplicities and their evolution with centrality at this energy (as seen in Figs. |T] and 
g), and the p^-spectrum in pp collisions at SPS and in pp collisions at SppS at CERN 
and TeVatron, and the increase of (pt) with energy and multiplicity (see [l^); we have 
also checked that this is neither an effect of pdf's or of their nuclear modifications, 
nor of rescattering, whose influence on the py-spectrum is tiny, see [|l6l and Fig. 
in fact, from the studies in |TB[ it can be concluded that the transverse momentum 
enhancement in collisions between nuclei compared to those between nucleons is due 
in the model both to the hard contribution which becomes more important with an 
increasing number of elementary collisions, and, above all, to the transverse momentum 
broadening of the partons at the ends of the strings introduced in the model and 
responsible of the increase of (pt) with increasing multiplicity, while string fusion has 
a very small effect. It is also remarkable that the discrepancy with the experimental 
data appears in a model like ours, which for the collisions studied at RHIC produces 
only 1 % of hard elementary collisions, and in a model like that of in which 
most of particle production at RHIC energies comes from the hard contribution^. So 
it really looks like an effect which diminishes the number of high pt partons, leading 



them to the low pt region. Jet quenching |3^ seems a good candidate to explain 
this experimental finding, but it should be taken into account that it also leads to the 
appearance of more particles at low px and rj; thus, the simultaneous comparison of the 
^Possible differences in tfie pr-spectrum in nucleon-nucleon collisions between our model and those 
based on hard scatterings like HIJING js^ should become visible at LHC, where the results are not 
so tightly constrained by the existing experimental data at SPS, SppS and TeVatron: In our model 
the contribution from hard scatterings will be smaller and thus we expect less high-pr particles. 
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evolution of both multiplicities and transverse momentum distributions with centrality 
should be a crucial test for this mechanism^. One would think that the presence of 
saturation of low transverse momentum partons p7| , ^ would make the comparison 
with experimental data even worse: the low pt region of the spectrum, populated of 
poorly resolved partons, would be damped due to parton fusion and the spectrum 
become flatter than without saturation. Quite the same would occur in percolation 
of strings soft strings have a larger transverse dimension than hard partons and 
would fuse more easily, and fused strings with higher string tension would produce 
particles with higher than ordinary strings, so the mean px would increase with 
atomic size or centrality contrary to what data apparently showQ. 

Finally, in Table |l| model results for different particle ratios are shown and com- 
pared with published experimental data [|T^, ^ ^ |^, ^, |T^, [TT|]. For completeness, 
let us indicate the results in the model for the ratios A/A, K^/K^, p/tt^ 

and K-/71- at r/ ~ 0, for which we get 0.85|0.87|0.87, 0.60|0.92|0.88, 1.08|1.03|1.04, 
0.02|0.07|0.04 and 0.08|0.12|0.16 respectively without string fusion or rescattering|with 
string fusion I with string fusion and rescatteringQ. The results in the model have been 
obtained in the corresponding pseudorapidity regions, for AuAu collisions at 130 GeV 
per nucleon in the center of mass with a centrality of 10 % and for particles with 
Pt > 0.2 GeV/c. Each experiment applies different centrality and kinematical cuts 
for the different ratios, but a common conclusion of all of them is that ratios are very 
weakly dependent on centrality of the collision and px of the particles, so this should 
not seriously affect the comparison. From these results it can be seen that the model 
overestimates antibaryon production, a feature already present at SPS, see [|1^, but 
string fusion is needed to increase the strangeness and antibaryon yield, which is badly 
underestimated, see the comparison with SPS data in M], if this mechanism is not 



^In 1^ the evolution of (j5}/(7r^) versus pr with centrahty is proposed as a test of jet quenching; the 

increase of this ratio with increasing pT observed by PHENIX Q is reproduced with a soft exponential 

component proportional to the number of participants plus a quenched perturbative distribution 

proportional to the number of binary collisions. In our model, the corresponding increase due to the 

soft part would be stronger than in due to string fusion and to the fact that this component is, 

in our case, proportional to the number of both wounded nucleons and binary collisions. 

recent analysis shows that nevertheless it is possible to simultaneously explain the evolution 

with centrality of both multiplicity distributions and transverse momentum spectra in a very crude 

realization of the percolating string approach. 

'^These results can be compared with preliminary, not yet published results: 0.73±0.03, 0.82±0.08, 

1.12±0.01±0.06, 0.08 and 0.15 respectively, presented by STAR at QM2001 fl. 
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included (in the ratios at central rapidities and due to the lack of stopping at RHIC 
energies, see below, and to the fact that string fusion creates on average the same 
amount of baryons and antibaryons, this feature is mainly visible in those involving 
multistrange baryons or in p/vr"). This discrepancy is less pronounced for S's than 
for A's, and for A's than for nucleons, and is more pronounced in the central region 
of (pseudo) rapidity. As stated in the brief model description, our rescattering model 
is simplistic, and cannot be expected to produce correct quantitative results, only the 
trend which it shows should be considered. So all that we can conclude is that for 
the ratios at RHIC, similar problems appear than those already present at SPS0. As a 
last comment, a preliminary, non-corrected for hyperon decay, measurement of the p-p 
yield at midpseudorapidity by BRAHMS 0, gives 8-t-IO for a centrality of 6 % (a value 
4-^6 has been extracted |Q from preliminary STAR data for the same centrality), 
while in our model we get a lower value ~ 2; this may suggest that the problem in the 
p/p ratio lies not only in a p excess, but also in some lack of stopping in the model. 

In conclusion, we have compared the results of the SFM with some of the first RHIC 
data. At RHIC, charged multiplicities in the central region for central collisions and 
their evolution with centrality are successfully reproduced, suggesting the presence of 
some mechanism, like string fusion, which moderates the increase of multiplicities with 
increasing centrality; On the other hand and in view of the SPS data, it is difficult to 
obtain clear conclusions from the behavior of multiplicities in the transition from SPS 
to RHIC. Results on particle ratios show, when compared to experimental data, similar 
problems of antibaryon excess previously found at SPS, and are probably related to 
the oversimplification of the model of rescattering and to problems with data at SPS, 
see [|16[. Finally, in the SFM the p^-spectrum at RHIC is fiatter than in data and this 
problem gets worse with increasing centrality, a feature which also appears in other 
models ^ in which the contribution of hard elementary collisions is much larger 
than in ours. At first sight, it looks improbable that parton saturation or percolation 
of strings could improve the comparison with the p^-distributions (but see ^1|]). So, 
from our point of view these data are most striking and, if confirmed, maybe a good 
candidate for a signature of non-conventional physics appearing in heavy ion collisions 
at RHIC. Although the results of the model on features which should depend strongly 
^Apparently, the antibaryon-to-baryon ratios measured at RHIC favor ^ a coalescence model 
||l|, see for a comparison of our results at full RHIC energy and those coming from other 
models. 
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on the evolution of the system (particle ratios and pT-spectrum if jet quenching is 
present) cannot be considered satisfactory, the agreement with multiplicities and their 
evolution with centrality, which are usually assumed not to vary too much during 



evolution gives us some confidence in the ability of the model to describe the 

initial condition, to be used for further evolution, in a collision between heavy ions at 
high energies. 
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List of figures 

1. Results of the model for the pseudorapidity distribution of charged particles for 
central (5 %) PbPb collisions at 17.3 GeV per nucleon in the center of mass (dashed- 
dotted line), and central (6 %) AuAu collisions at 56 (dotted line), 130 (dashed line) 
and 200 (solid line) GeV per nucleon in the center of mass, compared with experimental 
data at SPS from NA49 |2|, g] (black square) and WA98 |2§ (black, upward pointing 



triangle), and at RHIC from PHOBOS ||^, (black, downward pointing triangle for 
56 GeV, open circle for 130 GeV and black circle for 200 GeV), BRAHMS [g (open 
square) and PHENIX [Q (open triangle). 

2. Pseudorapidity density of charged particles at = divided by one half the number 
of participant nucleons, versus the number of participant nucleons, in PbPb collisions 
at 17.3 GeV per nucleon in the center of mass (multiplied by 1/2, lower curves and 
symbols) and in AuAu collisions at 130 GeV per nucleon in the center of mass (upper 
curves and symbols); also the experimental number for pp collisions at 130 GeV per 
nucleon is given [^, filled square. Experimental data are from PHENIX [Q (filled 
triangles), PHOBOS [0 (open triangle), WA98 H (filled circles) and NA49 |2|, § 
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(open circle). Curves are results of the model for the 75 % more central events, without 
fusion or rescattering (dotted lines), with fusion (dashed lines) and with fusion and 
rescattering (solid lines). 

3. The same as in Fig. 0, but for PbPb collisions at 5.5 TeV per nucleon in the center 
of mass. 

4. Transverse momentum spectrum (l/(27rpr) dN / {d7]dpT)\r]=o versus pt) of charged 
particles at ?7 = in AuAu collisions at 130 GeV per nucleon in the center of mass, for 
central collisions (5 %, solid and dashed lines and filled circles) and for minimum bias 
collisions (92 %, multiplied by 0.01, dotted and dashed-dotted lines and open circles). 
Data are from PHENIX 0] ; solid and dotted lines are results of the model with string 
fusion, dashed and dashed-dotted lines with string fusion and rescattering. 



List of tables 

1. Different particle ratios in central (10 %) AuAu collisions at 130 GeV per nu- 
cleon in the center of mass in the model without string fusion or rescattering (NF), 
with string fusion (F) and with string fusion and rescattering (FR) for particles with 
Pt > 0.2 GeV/c, compared with experimental data [|1^, H, 0, |^, |10|, |1^]. For the 



centrality criteria and kinematical cuts in the different experiments and ratios, see the 
experimental references and comments in the text. 
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Tables: 



Ratio 


NF 


F 


FR 


BRAHMS 


PHENLX 


PHOBOS 


STAR 


P/P 


0.81 


0.85 


0.80 


0.64±0.04 


0.64±0.01 


0.60±0.04 


0.65±0.01 


(r/~0) 








±0.06 {y ~ 0) 


±0.07 


±0.06 


±0.07 


p/p 


0.38 


0.50 


0.38 


0.41±0.04 








(l/~2) 








±0.06 








K-/K+ 


0.92 


0.97 


0.96 






0.91±0.07 




(77-0) 












±0.06 






1.02 


1.02 


1.01 






1.00±0.01 




(77-0) 












±0.02 





Table 1: 
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